A theoretical analysis of small Pt particles on rutile TiO2(110) surfaces 
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The adsorption profiles and electronic structures of Pt n (n=l-4) clusters on stoichiometric, 
reduced and reconstructed rutile TiO2(110) surfaces were systematically studied using on site d-d 
Coulomb interaction corrected hybrid density functional theory calculations. The atomic structure 
of small Pt cluster adsorbates mainly depend on the stoichiometry of the corresponding titania 
support. The cluster shapes on the bulk terminated ideal surface look like their gas phase low 
energy structures. However, for instance, they get significantly distorted on the reduced surfaces with 
increasing oxygen vacancies. On non-stoichiometric surfaces, Pt-Ti coordination becomes dominant 
in the determination of the adsorption geometries. The electronic structure of Pt„/TiO2(110) 
systems can not be correctly described by pure DFT methods, particularly for non-stoichiometric 
cases, due to the inappropriate treatment of the correlation for d electrons. We performed DFT+U 
calculations to give a reasonable description of the reconstructed rutile (110) surface. Pt clusters 
induce local surface relaxations that influence band edges of titania support, and bring a number 
of band-gap states depending on the cluster size. Significant band gap narrowing occurs upon Pt„- 
surface interaction due to adsorbate driven states on the bulk terminated and reduced surfaces. 
On the other hand, they give rise to a band gap widening associated to partial reoxidation of the 
reconstructed surface. No metallization arises even for Pt4 on rutile. 

PACS numbers: 68.43. Be, 68.43.Fg 



I. INTRODUCTION 

Transition metal oxides are important because of their 
wide range of technological applications. Naturally oc- 
curring rutile polymorph of titanium dioxide is a generic 
material because of its abundance, non-toxicity, and sta- 
bility under atmospheric conditions. The (110) termina- 
tion of rutile Ti02 is energetically the most stable sur- 
face among the other low index facets^ For these rea- 
sons, rutile TiO2(110) structure is a prototypical mate- 
rial to understand the catalysis on more complex oxide 
surfaces Its reducibility (either by oxygen vacancy for- 
mation or alkali metal adsorption) raises a great interest 
for fundamental study of photo- and heterogeneous catal- 
ysis^—, functional ultrathin films^, and dielectrics^^ 
These applications rely mainly on the surface electronic 
properties of titania. 

In Ti02 bulk structure, Ti atoms are sixfold and oxy- 
gens are threefold coordinated. Rutile (110) surface has 
these types of atoms, such as basal atoms Bl and Ti6c, 
which exhibit bulk like bonding characteristics as de- 
picted in Fig. [TJ The (110) termination of the bulk lat- 
tice breaks the Ti-0 bonds which lie normal to the sur- 
face plane resulting in fivefold Ti and twofold O atoms 
on the surface. Those are denoted as Ti5c and 01 in 
Fig. [JJ respectively. The undercoordinated bridging oxy- 
gens, 01, are exposed on the surface. They form oxygen 
rows along [001] direction. The atomic disposition of this 
stoichiometric long range ordered (lxl) phase has been 
well established by experimental methods^iii - — and by 
ah initio calculations *£r— 

By means of ion bombardment or thermal anneal- 
ing, the surface can be reduced through oxygen removal. 



Its oxygen vacancy induced reducibility enhances hetero- 
and photocatalytic activity by giving rise to rich surface 
chemistry. In this sense, rutile (110) system becomes an 
excellent model substrate that shows interesting proper- 
ties of more complex metal oxides^ - — 
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FIG. 1. (color online) Optimized atomic structure of bulk 
terminated rutile TiO2(110) surface. Actual computational 
cell has 15 atomic layers also referred as 5 trilayers. The model 
shows only 3 out of 5 trilayers. Black (red) and white balls 
represent O and Ti atoms, respectively. Dashed rectangular 
area denotes the (lxl) surface unit cell. The bond labels are 
indicated as capital letters that run from A to K. 
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Although DFT is successful in describing ground state 
atomic and electronic structure of many semiconductors, 
there is still a debate on the choice of method, pure or hy- 
brid, employed within the framework of DFT to improve 
the accuracy to estimate the material properties. Partic- 
ularly, some of DFT predicted properties of reconstructed 
and reduced Ti02 surfaces are inconsistent with the ex- 
perimental datai 41 ! 42 This is due to inadequate descrip- 
tion of strongly correlated 3d electrons of Ti atoms4£~— 

Since bulk termination gives rise to a sharp discontinu- 
ity in the atomic bonds, rutile (110) is known occasionally 
to undergo surface morphological reconstructions 
This phase transition is believed to be assisted by anneal- 
ing that reduces the surface and leads to a (1x2) recon- 
struction. The three dimensional determination of such 
atomic rearrangements on reconstructed rutile (110) sur- 
face is difficult by experimental techniques.— Onishi and 
Iwasawa proposed the Ti203 added row model for the 
reconstructed rutile (110) surface.— There are also some 
controversial experimental results. For example, Park et 
aL— suggested a new model where Ti are on interstitial 
sites forming a relatively more O deficient surface stoi- 
chiomctry. This model recently confirmed experimentally 
by Shibata et alM- Although it is a matter of debate, the 
added row model of Onishi and Iwasawa remains to be 
widely assigned (lx 2) rutile (110) reconstruction^ 7 -"— i5S 

The steps and point defects such as oxygen vacancies 
and (sub)surface impurities shown in scanning tunnel- 
ing microscopy (STM) images play an important role in 
electronic properties of TiC>2 surfaces 4 ^. Understanding 
of the fundamental aspects of defects and impurities on 
titania surfaces are important for technological applica- 
tions. In addition, the catalytic activity of Ti02 surfaces 
can be promoted by transition metal (such as Au or Pt) 
doping^ or adsorptio n 57 ' 58 . Therefore, the interaction 
of surface defects with functional adsorbates is crucial 
for both fundamental research and practical applications. 
For example, Pillay and Hwang 5 ^ studied the adsorption 
geometries of small Au, Ag, and Cu clusters on rutile 
(110) surface. Recently, Gong et al£^ investigated the 
atomic structures of small Au and Pt clusters supported 
on defect and regular sites of anatase TiO2(101) surface. 

In this work, our primary aim is to study the ad- 
sorption sites, geometries and resulting electronic struc- 
tures of Pt„ (n=l-4) clusters on the stoichiometric, 
reduced (by bridging oxygen Bl removal) and recon- 
structed (added row model 4 ^) rutile TiO2(110) surfaces 
using ab initio calculations. In order to gain insight 
into the growth pattern of small Pt particles on these 
surfaces, we chose a single Pt adsorbate as the starting 
point. Then, we studied the adsorption profiles of Pt„ by 
adding extra Pt atoms at different probable adsorption 
sites to the previously optimized Pt/TiO2(110) system. 
Alternatively, we also put Pt particles, whose geometries 
were optimized in their gas phase, on probable adsorp- 
tion sites. By comparing these two distinct approaches 
we decided on the lowest energy adsorption structures of 
small Pt particles on the surfaces considered. 



II. COMPUTATIONAL DETAILS 

Total energy density functional theory (DFT) calcula- 
tions were carried out using the Vienna ab-initio simu- 
lation package (VASP)i^ 1 - Nonlocal exchange-correlation 
energies were treated with Perdew-Burke-Ernzerhof 
(PBE)£2 functional based on the generalized gradient 
approximation (GGA). We used projector-augmented 
waves (PAW) method^ 3 -^ to describe the ionic cores and 
valence electrons with an energy cutoff value of 400 eV 
for the plane wave expansion. 

Pure DFT methods systematically fail in describing 
strongly correlated 3d electrons localized on Ti atoms, 
particularly in the case of reduced rutile (110) sur- 
face42~— This comes from the inherent shortcoming of 
the DFT due to the lack of proper self-energy cancellation 
between the Hartree and exchange terms as in Hartree- 
Fock theory. DFT+U approach is an alternative attempt 
to compensate this localization deficiency by introducing 
additional Hubbard U term for electron on-site repul- 
sion. Therefore, in order to get a sound description of 
Pt n adsorption on reduced and reconstructed TiO2(110) 
surfaces, we performed Hubbard U corrected DFT cal- 
culations using Dudarev's approach^ We supplemented 
PBE exchange-correlation functional with Dudarev U—J 
term acting on Ti 3d states, that we refer as GGA+U. 
This empirical term is adjusted to give reasonable agree- 
ment with the experimental data. Morgan and Watson 
fitted the U parameter to 4.2 eV in order to reproduce 
the experimentally observed position of oxygen vacancy 
driven Ti defect state in the band gap 4 s - In the same 
manner, our tests suggested to choose the parameters 
[7=4.5 eV and J=0 eV to obtain a reasonable electronic 
description of the added row model of Onishi and Iwa- 
sawa as well as the gap state at the reduced surface. 

Theoretical studies employing other hybrid methods 
such as incorporation of the exact Fock exchange energy 
(EXX) with various percentages, recently, revisited the 
stoichiometric rutile TiO2(110) surfaced Zhang et al. 
reproduced many of the experimental bulk and stoichio- 
metric surface properties of Ti02 by adjusting the mixing 
of HF exchange contribution. Further investigations are 
still needed to see if EXX approach can correctly describe 
the defect states. 

Computationally, there has been another debate about 
the slab thickness to correctly predict the surface ener- 
getics as well as the atomic structure. The surface energy 
is known to fluctuate with slab thickness ^ 19 ' 28 ' 31 " — For 
adsorption of small molecules on the bulk terminated sur- 
face, Thompson et al£^ suggested 4 layers of 0-Ti2 02- 
O units (also referred as trilaycrs) with the two bottom 
layers fixed to their bulk positions. Kowalski et alM- ad- 
ditionally suggested the saturation of the dangling bonds 
at the bottom surface with pseudohydrogens having nu- 
clear charges of +4/3 and +2/3, correspondingly. 

Our stoichiometric slab model consists of 15 atomic 
layers that correspond to 5 layers of 0-Ti202~0 units. 
This model has a mirror symmetry with respect to the 
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central atomic plane. The dashed rectangular region 
shown in Fig. [TJ represents p(l x 1) periodicity. The cor- 
responding surface unit cell is composed of a central 
Ti2C>2 planar arrangement with two oxygens symmet- 
rically cross bonded to the nearest neighbor in-plane Ti 
atom from above and from below. These units are shifted 
by half its length along [110] between adjacent layers as 
seen in Fig. [T] which shows three out of the five layers in 
our model. The slab was separated from its images along 
the surface normal by a vacuum region of -14 A. 

For geometry optimizations, the Brillouin zone integra- 
tions were carried out with 2x2x1 Monkhorst-Pack^ 
/c-point grid for unreconstructed (bulk terminated and 
reduced) and reconstructed cells with surface periodic- 
ities p(4x2) and p(4xl), respectively. The spin polar- 
ization was found to be negligibly small for the Pt„ on 
stoichiometric surface considered in this work. This is 
in agreement with the similar findings of Iddir et al. 
for single Pt adsorbate on TiO2(H0)££. We performed 
spin polarized GGA+U calculations for Pt„/TiC>2 sys- 
tems where magnetization is non-negligible, in particular, 
for reduced and reconstructed surface cases. A full geom- 
etry optimization was fulfilled with residual minimization 
method direct inversion in the iterative subspace (RMM- 
DIIS) scheme preconditioned by a few non-selfconsistent 
Davidson-block iterations as implemented in the code. 
We required a precision of 10~ 2 eV/A in the residual 
forces in every spatial component on all the atoms with- 
out fixing them to their bulk positions. 

We estimated the binding energies of Pt clusters by, 

Eptv = £pt n /Ti0 2 — E>Ti0 2 ~ -Ept„ j 

where E Ptn / Ti Q 2 , i?Tio 2 an( i ^Pt„ are the total energies of 
the Pt„/TiC>2 combined system, the bare TiC>2 slab and 
the corresponding Pt„ cluster, respectively. The average 
adsorption energy per atom, Ep t ^ , can be calculated by 
dividing Ep t by the cluster size n. 



III. RESULTS & DISCUSSION 

We made extensive tests for the choice of the slab 
model by changing the slab thickness, the number of 
fixed atomic layers within, and the size of the vacuum 
spacing, provided that the (110) facet possesses 4x2 cell 
symmetry to avoid any interaction between the periodic 
images of Pt cluster adsorbates. The results suggest that 
no atom, particularly near the surface, should be fixed 
to its bulk position in order to avoid stress driven gap 
states and that the slab model has to be at least 5 tri- 
layers thick (15 atomic layers) to ensure a bulk-like cen- 
tral part. Similarly, for instance, Kiejna et al£^ found 
that surface energetics are very sensitive to lattice re- 
laxations by considering the effect of various slab model 
thicknesses. Other hybrid DFT approaches with different 
exchange and correlation treatments are known to have 
a small effect on the surface atomic configurations 3 ^. 



TABLE I. The comparison of computational and experi- 
mental bond lengths for the bulk terminated stoichiometric 
TiO2(110) surface. All measurements are in angstroms. Bond 
labeling follows similar to those of Thompson et al. as shown 
in Fig. [TJ 

Bond Experimental Theoretical 
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Ref.— 


Ref.— 


Ref.^ Ref.^. 


GGA_^ 


GGA+U 1 


A 


1.71±0.07 


1.85 


1.80 


1.84 


1.85 


1.88 


B 


2.15+0.09 


2.15 


2.04 


2.04 


2.03 


2.03 


C 


1.99+0.09 


2.08 


2.09 


2.11 


2.10 


2.09 


D 


1.84±0.05 


1.90 


1.95 


1.92 


1.95 


1.95 


E 


1.84+0.13 


1.79 


1.85 


1.83 


1.84 


1.90 


F 


1.97+0.12 


1.90 


1.90 


1.89 


1.90 


1.95 


G 


1.99+0.05 


2.00 


1.97 


1.98 


1.98 


1.97 


H 


2.18±0.11 


2.11 


2.11 


2.13 


2.12 


2.12 


I 


2.00±0.08 


2.01 


2.02 


2.02 


2.02 


2.01 


J 


1.92±0.06 


1.92 


1.97 


1.96 


1.96 


1.96 


K 


1.94±0.06 


1.89 


1.90 


1.91 


1.92 


1.95 



a values for 5 layer slab model 

b results of present calculation with PBE xc functional 

Indeed, a comparison of our calculated bond lengths 
for the stoichiometric rutile TiO^llO) surface shows a 
good agreement with the available experimenta l 1 5 ' 1 7 and 
theoretica l 25 ' 32 results as presented in Table HI Moreover, 
the incorporation of supplemental Coulomb repulsion to 
correct the correlation energy by inclusion of Dudarev 
U = 4.5 eV term for Ti 3d electrons did not distort 
the atomic positions, considerably. In addition to these 
structural properties, we also calculated the binding en- 
ergies (BE) of Pt to be 2.32 eV and 3.32 eV on the sto- 
ichiometric and partially reduced (with oxygen vacancy 
concentration of 9=1/8) surfaces, respectively. Our val- 
ues are presented in Table [TT] and are only slightly larger 
than the previous theoretical results of Iddir et al£2. This 
can be well addressed to the difference in slab thicknesses 
and in surface areas of computation cells. In particu- 
lar, the oscillatory convergence behavior of the surface 
electronic properties depending on the number of trilay- 
ers, included in a slab calculation, was also confirmed 
by other studies! 31 ' 32 ' 34 ' 35 From computational point of 
view, a slab with large number of trilayers is expected to 
give very well converged results for Pt clusters on 4x2 
cell. However, independent of slab thickness, the elec- 
tronic description of defect states such as the oxygen va- 
cancies is problematic with pure DFT methods^ 3 - - — In 
the same line, our tests showed that standard DFT fails 
in describing the electronic structure of the added row 
model of reconstructed surface, as well, by giving surface 
Ti 3c? state in the conduction band (CB). This study, par- 
ticularly, focuses on a reasonable electronic description of 
small Pt clusters on rutile (110) surface where spin polar- 
ization is non-negligible. Therefore, our GGA+U values 
for this slab model would not only give very well con- 
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verged bond lengths, as presented in Table |H but also 
yield the energetics reasonably accurate for a Pt clus- 
ter adsorption system on specifically non-stoichiometric 
rutile (ffO) surfaces. 

In order to study the adsorption characteristics of Pt„ 
(n = f - 4) particles on rutile TiC^llO) surfaces we first 
considered these in the gas phase and performed non- 
collinear spin polarized DFT calculations. A recent work 
stresses the necessity for inclusion of spin-orbit interac- 
tions in a relevant calculation to understand the physics 
of Pt clusters^ 7 -. However, the ground state geometries 
being a dimer for Pt2, a triangle for Pt3 and a bent (off- 
planar) rhombus in the case of Pt4 seem to be less af- 
fected by spin-orbit coupling. 

For Pt dimer, our non-spin-orbit calculations resulted 
in triplet electronic ground state with a binding energy 
(BE) of 1.819 cV/atom and gave the bond length to 
be 2.331 A. The calculated dimer length shows an ex- 
cellent agreement with the experimental value of 2.333 
A— — , while our result for the BE is slightly higher 
than the experimental value of 1.570 eV/atom. When 
we included self-consistent noncollincar spin-orbit cou- 
pling that resulted in the same magnetic ground state, 
the BE decreases to 1.665 eV/atom whereas the dimer 
bond length increases to 2.382 A. These values are sen- 
sitive not only to the inclusion of spin-orbit effects but 
also to the exchange-correlation scheme employed. For 
instance, by performing GGA-PW91 calculations, Huda 
et al£l determined different spin multiplicities for the 
ground state of Pt3 cluster depending on the spin-orbit 
coupling. However, our GGA-PBE calculations gave sin- 
glet ground states for Pt3 having average binding ener- 
gies of 2.184 and 2.376 cV/atom with and without LS 
coupling, respectively. Moreover, Huda et al. reported 
that spin-orbit interaction drives the geometry from an 
equilateral to an isosceles triangle. Pt3 with GGA-PBE, 
on the other hand, develops equilateral coordination with 
bond lengths of 2.49 A while the inclusion of spin-orbit 
interaction extends all bonds slightly and equally to 2.50 
A. Our GGA-PBE binding energy value agrees well with 
a previous theoretical result.— In addition, incorpora- 
tion of spin-orbit coupling leads to an excellent agreement 
with an old experimental estimate of 2.18 eV/atomi 7 - 2 - 

Pt4 has a quintet electronic ground state (GS) with 
GGA-PBE, leading to an off-planar rhombus geometry 
that has a side bond length of 2.51 A and a bending 
angle of 23.3°. We calculated atomic BEs to be 2.515 
eV/atom with, and 2.686 eV/atom without the LS cou- 
pling. These values are comparable with the previous re- 
sults i^liZ 3 . However, spin-orbit interaction included in the 
GGA-PBE calculation did not make the geometry per- 
fectly planar as reported for GGA-PW91.— It slightly 
extends the bond lengths to 2.52 A by decreasing the 
bending angle to 13.1° and gives the same spin multi- 
plicity for the GS. Consequently, for all isolated small Pt 
particles considered in this study, inclusion of LS coupling 
in GGA-PBE exchange-correlation scheme does not de- 
form the GS cluster shapes into new geometries. 



A. Pt n on stoichiometric rutile TiO2(110) surface 

Pure DFT with GGA-PBE functional gives a band gap 
of 1.48 eV for the stoichiometric rutile TiO2(110) sur- 
face. Supplemental Hubbard U on-site repulsion, acting 
on the Ti 3d electrons, heals this band gap to 2.01 eV by 
correcting the underestimation in the correlation energy. 
Although the Dudarev U parameter can be adjusted to 
a larger value to get a better agreement with the exper- 
imental value of ~3 eV— , this would significantly distort 
the atomic structure. Therefore, without losing the con- 
sistency of calculated atomic positions and bond lengths 
with the experimental values (as presented in Table [TJ) , 
we set it to U = 4.5, in consistency with the previous 
theoretical studies^— , to get the Ti defect state in the 
band gap of the reduced surface and to obtain a reason- 
able electronic description of the added row model of the 
reconstructed surface. In both reduced and reconstructed 
cases, pure DFT methods incorrectly pin the defect states 
inside and to the lower part of the conduction band (CB) 
which indicates a tendency for strong derealization. 

GGA+U calculations were performed using 4x2 cell 
in order to avoid any charge transfer between the peri- 
odic images of Pt clusters on the bulk terminated (110) 
surface. Moreover, in order to obtain the energetics 
correctly, we first investigated the magnitude of spin- 
polarization for small Pt clusters on the surface. The 
results gave negligibly small spin multiplicities for the 
ground states of TiO2(110) supported Pt„ clusters and 
for the bare surface, as well. In particular, similar find- 
ings were reported by Iddir et al. for Pt/TiO2(110) sys- 
tem^ Therefore, adsorption profiles [Fig. [5] and elec- 
tronic structures [Fig. [3] were calculated without spin 
polarization for the case of stoichiometric surface. 

Relaxed atomic positions for defect free surface are pre- 
sented in Tabic |H Our GGA and GGA+U values do not 
considerably differ from each other and are in good agree- 
ment with previous theoretical calculation o 22 ' 25 ' 29 ' 31 " — 
as well as the experimental results! 15 ' 17 The only discrep- 
ancy which exists for all DFT studies, is seen in Ti6c-Bl 
bond (B in Fig. [l}. In order to inspect its relevance to 
the number of layers, we also performed relaxation cal- 
culations with 7 and 8 trilayer cells, that reproduced the 
same value. Therefore, it can be well addressed to GGA- 
PBE exchange-correlation functional that underestimcs 
Ti6c-Bl in-plane bond length by 0.12 A. In addition, 
experiments reported slightly different results for bond 
length between Ti6c and the nearest neighbor subsurface 
oxygen (C in Fig. [T]). For this particular bond distance, 
our GGA and GGA+U values are consistent with the re- 
sult of Charlton et alr^ more than the recent result of 
Lindsay et al^ Agreement with the experimental data 
gets better for deeper subsurface bond lengths. 

The density of states (DOS) structure for the bare sur- 
face is presented in the first panel of Figj3] The Fermi 
energy is set to be the zero of energy scale and located 
just above the valence band maximum (VBM). The up- 
per part of the VB shows dominant O 2p character with a 
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FIG. 2. (color online) Minimum energy structures of Pt„ 
(n = 2-4) adsorbed on stoichiometric rutile TiO2(110) sur- 
face. O, Ti, and Pt atoms are denoted by black (red), white 
small balls and gray big balls, respectively. Pt dimer, trimer, 
and tetramer on this surface are depicted in (a), (b), and (c), 
respectively. Measurements are shown, on the corresponding 
bonds, all in angstroms. The top view of the small Pt parti- 
cles are presented in the insets to provide better description 
of the adsorption sites. 



non-negligible Ti 3d contribution. Top of the VB mainly 
composed of 01 states giving a peak ^0.4 eV below the 
VBM. The contribution from the basal oxygens, Bl, lies 
relatively lower in energy inside the VB. Our site pro- 
jected DOS analysis for the VB agrees well with previous 
studies, (e.g. hybrid PBEO result of Labat et al^ and 
pure GGA result of Sano et al*&). The lower part of the 
conduction band (CB) is dominated by Ti 3d states. The 
bottom of CB shows a hybrid Ti5c-Bl character with a 
large Ti contribution. Although, the composition of the 
CB have some similarities, we have considerable disagree- 
ment in the shape of the CB edge with other GGA and 
PBEO result o 29 ' 35 . This is due to on-site U repulsion act- 
ing on the Ti 3d electrons that causes charge localization 
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FIG. 3. (color online) Calculated projected and total density 
of states (DOS) structures for Pt„ (n= 1-4) on stoichiomet- 
ric rutile TiO2(110)-4x2 surface. Upper three panels corre- 
spond to the atomic structures presented in Fig. [2] 



around the atomic Wigner-Seitz radii. Dudarev's U sup- 
plements correlation energy of d electrons that elevates 
the corresponding band offsets to higher energies giving 
a gap of 2.01 eV. Moreover, it induces differences in the 
dispersion of Ti 3d states, particularly, near the conduc- 
tion band minimum (CBM). 

Binding of a single Pt atom to the stoichiometric sur- 
face was investigated for all possible adsorption sites. We 
found the lowest energy position of Pt at the hollow site 
over the Ti5c, leaned toward the nearest bridging oxy- 
gen, 01. This adsorption position was also reported pre- 
viously by GGA-PBE calculations^. The calculated BE 
is 2.32 eV at the hollow site where Pt-0 and Pt-Ti bond 
lengths become 1.96 and 2.49 A, respectively. Pt atom 
pulls Ti5c slightly up out of the surface plane by 0.15 
A, and also draws 01 off the bridging oxygen row so 
that the bond length between this oxygen and Ti6c ex- 
tends to 2.04 A from its bare surface value of 1.88 A [in 
Table H]. These values are slightly larger than the GGA- 
PBE results due to additional on-site Coulomb repulsion 
between Ti 3d electrons with Dudarcv U — 4.5 eV. For 
instance, 01-Ti6c bond length increases by 0.03 A from 
its GGA-PBE value of 1.85 A. 

The DOS for Pt/TiO 2 (H0) presented in Fig. [3] indi- 
cates that single Pt at the hollow site interacts with Ti5c 
that appears as a Pt-Ti5c hybridization peak at the bot- 
tom of the CB. The non-bonding excess charge around 
Pt brings four occupied flat gap states within a range 
of 0.9 eV above the VBM. In addition to the band gap 
narrowing, these impurity states increases visible range 
transition probability upon a vertical excitation. This 
will lead to a metal-to-substrate charge transfer, useful 
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for solar cell applications. 

Energetically preferential adsorption geometry of Pt 
dimer is shown in Fig. [5^,. This atomic structure happens 
when a second Pt is attached to the Pt/Ti02 system as 
to form a dimer. While one of the Pt atoms is already at 
the hollow site the other one finds its minimum energy 
position above basal oxygen, Bl. It also weakly pulls 
the nearest neighbor 01 slightly distorting it out of the 
oxygen row at a separation of 2.39 A. Alternatively, we 
placed Pt dimer as a whole on the surface at possible 
adsorption sites and found the same geometry between 
the oxygen rows as the minimum energy structure. Pt 
Pt bond distance becomes 2.57 A that is meaningfully 
larger than the isolated Pt dimer length by 0.24 A. 

Adsorption of Pt trimer breaks its gas phase equilateral 
symmetry and causes local deformations at the surface 
plane and in the second subsurface layer as seen in Fig [2b. 
This geometry is obtained through Pt2/Ti02 system by 
an additional Pt atom as to form a triangular clustering 
so that it coordinates with the nearest neighbor 01 at 
the second oxygen row. Pt3 at different adsorption sites 
is energetically not preferable. The two Ti5c under the 
cluster are lifted up by ^0.23 A and Bl atoms, coordi- 
nated to these Ti5c, are pushed down by 0.5 A from their 
relaxed clean surface positions. The longest Pt-Pt bond 
length occurs as a result of the pulling of the bridging 
oxygens from the two sides so that it slightly stretches 
from 2.50 A to 2.56 A. Meanwhile, the shortest Pt-Pt 
bond is 2.46 A between the Pt's at the hollow and the 
Ti5c sites. 

Electronically Pt2 and Pt3 on the stoichiometric sur- 
face yield very similar DOS structures except the number 
and the band energies of cluster-driven impurity states 
below the Fermi energy. Pt dimer brings six occupied 
states within the range of 1.30 eV above the VB and an 
empty flat going gap state 0.15 eV below the CB. There- 
fore, Pt2 on the surface causes a band gap narrowing of 
1.45 eV with respect to that of the clean surface. The 
interaction between the additional Pt at Pt3 and the sub- 
strate occurs essentially with the nearest neighbor 01 on 
the second oxygen row in Fig. It brings two more oc- 
cupied flat states, one of them at the VBM, and causes 
a shift in band energies of the gap states. In the case of 
Pt3, metal driven empty gap state lies almost at the same 
position with that of Pt 2 /TiO2(110). Resulting band gap 
narrowing is calculated to be 1.50 eV. 

Pt-Ol coordination number increases to three as the 
cluster size adds up to Pt4. Although it gives the 
strongest total binding of 2.98 eV (in Table HI)) . BE per 
Pt of 0.75 eV proves to be the lowest among the other ad- 
sorbates. The reason for this is that one of the Pt atoms 
lying above the other three makes no contact with the 
surface. Pt4 particle on the surface occurs to be slightly 
distorted bent-rhombus similar to its isolated GS struc- 
ture (in Fig. Ht). Electronically, Pt 4 /TiO 2 (110) system 
features the lowest energy band gap by filling the upper 
lying flat gap state that is 0.13 eV below the CB. The 
position of this state is almost the same for the multi- 



platinum adsorbates and is empty in the cases of Pt 2 
and Pt3 structures. Altogether, Pt 4 brings a total of ten 
occupied gap states. It couples to the surface by giv- 
ing significant contribution to the upper part of the VB 
essentially by three Pt-Ol hybrid bondings and at the 
bottom of the CB by three Pt-Ti5c coordinations. 

For all cases, adsorption of small Pt particles causes 
considerable local distortions that are mediated by the 
neighboring atoms to the second trilayer underneath and 
are proportional to the cluster size. Although Pt„ ad- 
sorbates induce surface deformations, clusters themselves 
seem to be less affected by forming structures similar to 
their gas phase low energy geometries on stoichiomet- 
ric rutile (110) surface. This indicates that the metal- 
metal coordination within the clusters is stronger than 
the cluster-substrate interaction. 



B. Pt n on reduced rutile TiO2(110) surface 

Reduced surface is constructed using the relaxed stoi- 
chiometric 4x2 cell with 5 trilayers. We removed a single 
bridging oxygen, 01, from both surfaces to avoid forma- 
tion of an unreal dipole across the slab. Isolated oxygen 
vacancies are experimentally observed to produce defect 
states in the band gap 0.7-0.9 eV below the CBMr22. 
Pure DFT methods tend to give too delocalized solution 
where excess charge density occupies the bottom of the 
CB. Morgan et alM suggested that a spin polarized GGA 
calculation, corrected by Hubbard U with U ^ 4.2 eV 
reproduces the experimentally observed gap state. Our 
choice of U = 4.5 eV falls within this range and gives the 
defect state in the gap ~1.2 eV below the CB as shown 
in the lower panel of Fig. [5j Spin multiplicity at the GS 
happens to be 4 for the 4x2 cell with two reduced sur- 
faces. The position of the defect states is sensitive to 
atomic relaxations that depends on the choice of the Du- 
darev U parameter. Calzado et al£^ also showed that 
gap states move closer to the CB edge as the oxygen va- 
cancy concentration increases in agreement with the UPS 
observation :) 36 ' 74 . 

Single Pt atom finds its minimum energy position at 
the bridging oxygen vacancy site with equal Pt-Ti6c 
bond lengths of 2.41 A each. Pt oxidizes the surface 
such that the vacancy induced topographical deforma- 
tions recover. The combined system looks like the 4x2 
stoichiometric surface cell except that one of the bridging 
oxygens replaced by a Pt atom which stays 0.7 A above 
the oxygen row level. Therefore, Pt adsorption at the 
defect site gives a strong binding of 3.32 eV compared 
with that of the stoichiometric case. 

We also investigated Pt adsorption on fully reduced 
rutile(llO) surface. At high oxygen vacancy concentra- 
tions, an isolated Pt atom takes the place of the nearest 
neighbor threefold coordinated basal oxygen, Bl; an ef- 
fect also known as the strong metal support interaction 
(SMSI). Then, Bl loses its coordination with Ti5c and 
moves to the nearest vacancy site to form a bridge con- 




FIG. 4. (color online) Minimum energy structures of Pt„ 
(n = 2 - 4) adsorbed on partially reduced rutile Ti02 (110) sur- 
face. O, Ti, and Pt atoms are denoted by black (red), white 
small balls, and gray big balls, respectively. Pt dimer, trimer, 
and tetramer on this surface are depicted in (a), (b), and (c), 
respectively. Measurements are shown, on the correspond- 
ing bonds, all in angstroms. The top view of the small Pt 
particles are presented in the insets are provided for visual 
convenience. 



figuration so that the surface gets oxidized. Pt at Bl site 
causes considerable local geometric deformations while 
Bl becomes a twofold coordinated 01. 

The site projected DOS structure of the Pt/pr-Ti02 
system is shown in Fig. [5] The contribution of single Pt 
adsorbate is seen below the Fermi energy indicated as the 
dark shade (red). The ground state of a Pt atom is cPs 1 . 
Therefore, Pt makes two strong bonds with Ti6c at the 
vacancy site. The BE is 3.32 eV this is 1 eV larger than 
that on the stoichiometric surface. These bonding energy 
states fall within the valence band. Pt at the defect site 
has a significant influence on the local surface reconstruc- 
tion. This is clearly seen when one compares the DOS 
structures at the VBM, before and after Pt adsorption 
on pr-Ti02- Pt adsorbate also brings four occupied gap 
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FIG. 5. (color online) Calculated projected and total DOS 
structures for Pt n (n = 1 - 4) on partially reduced (pr) rutile 
TiC>2(110)-4x2 surface. Upper three panels correspond to the 
adsorption geometries in Fig. [4] 



states. Upper lying two of them have strong Pt contri- 
bution at 0.78 and 1.37 eV above the VBM. Only the 
third one below the Fermi energy shows a weak disper- 
sion. Others are all flat showing non-bonding character. 
Lowest lying flat band is 0.43 eV above the VBM. 

Adsorption geometries of Pt„ for n = 2 - 4 on oxygen 
defective surface are shown in Fig.@]whcrc pr-TiO2(110) 
stands for partially reduced rutile TiO2(110) surface as- 
sociated with oxygen vacancy concentration of 1/8. Pt2 
binds to the reduced surface at the vacancy site as shown 
in Fig. H^i. Each Pt atom relaxes atop Ti6c site at each 
side so that the dimer tends to align with the oxygen row. 
Pt2 at the defect site pushes Ti6c atoms down from their 
in plane positions causing considerable distortion on the 
Ti6c row. Pt-Pt bond length shortens to 2.40 A from 
its isolated value of 2.52 A as a result of strong binding. 
The BE per Pt atom was calculated to be 1.11 eV that is 
larger than that of Pt dimer on the stoichiometric surface 
due to the oxygen vacancy. 

We have traced a number of possible adsorption ge- 
ometries for Pt3 particle on the defective surface. Pta/pr- 
TiO2(H0) in Fig. [4]3 represents the minimum energy ge- 
ometry of the trimer whose initial structure is obtained 
by adding a Pt atom to the Pt dimer in Fig. This 
additional Pt gets adsorbed at Ti5c site and pulls it up 
by ^0.40 A from its in plane position. Resulting Pt-Ti5c 
bond length becomes 2.42 A. Due to imperfect alignment 
of the Pt dimer along the surface oxygen row this trian- 
gular cluster shape becomes slightly distorted compared 
with the equilateral triangle of gas phase Pt3 . 

Pt2 and Pt3 adsorbates give very similar DOS charac- 
teristics on the reduced surface. For both cases, eight im- 
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purity states fall within the band gap as shown in Fig. [5] 
One expects to get narrower band gaps as the number 
of Pt on the surface increases. Interestingly, Pt2/pr- 
TiO2(H0) system has the narrowest band gap of 0.05 
eV among the other cases. This corresponds to a band 
gap narrowing of 1.22 eV relative to the Fermi energy 
of the pr-TiO2(110). Although the band gap is under- 
estimated by DFT due to improper description of the 
exchange-correlation energy, band gap narrowing must 
be absolute and so, experimentally verifiable. This value 
has been calculated to be 1.16 eV for the Pt3 case. 

In Pt4 case, we examined initial adsorption structures 
including low energy gas phase rhombus and pyramid 
isomers on partially reduced surface. Relaxed cluster- 
surface combined system has been found to be the geom- 
etry as shown in Fig. HJ;. Symmetrical rhombus form can 
be obtained by adding a Pt atom to Pt3 / pr-TiC>2 system 
in Fig. |4Jd. However, after relaxation, such a structure 
is energetically less favorable by 0.35 eV/cell. Compara- 
bly, Pt4 in Fig. 2t, shows an increased coordination with 
Ol's. This structure can also be obtained by adding a Pt 
to Pta/pr-Ti02 so that Pt4 lies between the two bridging 
oxygen rows. All these small Pt particles cluster around 
the vacancy site. Their geometries are different, in par- 
ticular for Pt4, from those on the stoichiometric surface, 
due to the surface charge distribution around the defect 
site. Pt particles oxidize and therefore bind to defective 
surface stronger than to the defect free one. 

Pt4 adsorption on the reduced surface induces thir- 
teen flat-like occupied gap states due to the excess charge 
brought by the metal cluster, which is larger than that 
of the previous cases. Fermi energy occurs at just above 
the upmost lying flat state as seen in the top panel of 
Fig. [5] This corresponds to a narrowing of 1.06 eV rel- 
ative to the band gap of the reduced surface. Since all 
Pt clusters bind so as to center around the defect site, 
all cases have contributions to the DOS'es from these 
bonding Pt-Ti6c states that appear inside the VB with 
increasing Pt content. Indeed, Pt4/pr-TiO2(110) system 
has the largest partial DOS contribution inside the VB 
which essentially disperses similar to that of the Pt3 case. 



TABLE II. Calculated total binding energies (eV) of Pt clus- 
ters and corresponding average binding energies per Pt atom 
at TiO 2 (110)-4x2 surfaces. 





stoichiometric 


reduced^ 


reconstructed^ 


Pt 


2.32 / 2.32 


3.32 / 3.32 


3.94 / 3.94 


Pta 


1.74 / 0.87 


2.22 / 1.11 


4.20 / 2.10 


Pt 8 


2.66 / 0.89 


3.26 / 1.09 


4.76 / 1.59 


Pt 4 


2.98 / 0.75 


3.38 / 0.85 


4.54 / 1.14 



a with oxygen vacancy concentration 9=0.125 
b added-row model of Onishi and Iwasawa2£ 



C. Pt n on reconstructed rutile TiO2(110) surface 

The atomic resolution of the long range (1 x 2) 
phase of reconstructed rutile TiO 2 (110) surface is dif- 
ficult with experiments, numerous theoretical studies at- 
tempted to resolve the controversy regarding the sur- 
face morphology 4 ^—. Scanning tunneling microscopy 
(STM) and low-energy electron diffraction (LEED) ex- 
periments showed the existence of Ti2 03 added rows in 
agreement with the model proposed by Onishi and Iwa- 
saw a 46 ' 48 ' 51 " — . We constructed this model of (1 x 2) sur- 
face by forming additional Ti203 rows along [001] at the 
top and the bottom facets of the relaxed stoichiometric 
slab with 7 trilayers. This turns out to be a symmetrical 
cell which prevents a fictitious dipole formation across 
the slab. We studied Pt cluster adsorption on 4 x 2 cell 
that accommodates separations of ~ 12 A along [001] and 
- 13 A along [010] to avoid any charge transfer between 
the periodic images of the metal adsorbatcs. The relaxed 
geometries of Pt„/ar-TiO2(110) systems are depicted in 
Fig. [6] 

The electronic structure of added row model of the 
reconstructed surface was given by Kimura et al£^ and 
recently by Blanco-Rcy et aL— . By employing pure DFT 
methods they obtained DOS peaks resulting from Ti203 
row occupying the bottom of the CB. The origin of these 
DOS contributions are similar to Ti +3 states derived from 
oxygen vacancies on the reduced surface. Pure DFT is 
known to fail in describing these vacancy states predict- 
ing a metallic character42r— However, defect states ap- 
pear 0.7-0.9 eV below the Therefore, for (1x2) 
reconstruction, pure DFT must incorrectly predict Ti 3c? 
states coming from the Ti2 03 row to appear inside, at 
the bottom of, the CB. We tried to obtain a reasonable 
electronic description of the added row model of recon- 
structed, bare and Pt„ (n=l-4) adsorbed surfaces by 
employing Dudarev U corrected spin polarized DFT cal- 
culations. Our tests showed that calculated band struc- 
ture of Ti02 is sensitive to relaxation of the atomic posi- 
tions based on energy minimization of the computational 
cell. In this sense, supplementary U repulsion energy 
considerably changes the Ti-0 bond lengths, resulting 
surface geometry, and so, the band structure. 

Experimentally, Abad et al. observed the band gap 
state at 0.7 eV having Ti 3d character from the UPS He- 
I spectra for TiO2(110)-(lx2)Z&. In our calculations, the 
choice of U = 4.5 drives the Ti defect states from the CB 
edge down into the gap giving band gap of ~0.1 eV in 
the bottom panel of Fig. [7j Certainly, this value is un- 
derestimated due to the improper cancellation of the self 
interaction, a well known artifact of DFT. In addition to 
the appearance of defect states due to broken stoichiom- 
etry, reconstruction also changes the dispersion of states 
at the band edges. Our GGA+U calculations resulted in 
a ground state with spin multiplicity of 1.99 per (1x2) 
surface cell similar to the case of pr-TiO2(110). Two 
satellite DOS peaks appear in the band gap associated 
with the excess charge on the Ti2 03 row. These defect 



states strongly disperse over a width of ~1.3 cV in the 
gap and possess Ti 3d character. Fermi level occurs 2.16 
eV above the VBM. If we compare band gap narrowing of 
2.16 eV with experimental gap, we expect (1x2) phase 
of TiO2(110) to exhibit a photoemission at an energy 
consistent with the observation of Abad et a/.— 







FIG. 6. (color online) Minimum energy structures of recon- 
structed rutile TiO2(110) surfaces with Pt„ (n = 1-4) adsor- 
bates. O, Ti, and Pt atoms are denoted by black (red), white 
small balls, and gray big balls, respectively. Bond lengths are 
shown in angstroms. The top view of the small Pt particles 
are presented in the insets are provided for visual convenience. 
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FIG. 7. (color online) Calculated projected and total DOS 
structures for Pt„ [n = 1 - 4) on added row (ar) model of re- 
constructed rutile TiO2(110)-4x2 surface for the low energy 
adsorption models shown in Fig. [6] 



We considered various possible sites for Pt adsorption 
on the reconstructed (1x2 reduced) surface. The min- 
imum energy surface structure of single Pt adsorption 
is depicted in Fig. [(Jl An isolated Pt atom gets ad- 
sorbed at the oxygen site over the terrace in between the 
two Ti cations lying along [010] that are also labeled as 
Ti(iv)^. Pt pulls these Ti cations, and so, indirectly 
the oxygen underneath, forming a symmetrical tctrago- 
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nal group. The oxygen elevates up by 1.5 A from its re- 
laxed Ti 2 03 row position. Pt-Ti equilateral bond length 
becomes 2.44 A and Ti-0 bonds are 1.89 A belonging to 
the tetragon while the two Ti(iv)'s make four equidistant 
bonds with oxygens over the row reading 1.87 A each. Pt 
on reconstructed surface has the strongest binding energy 
of 3.94 eV/atom among the other system considered in 
this study. 

Electronically, adsorbate driven atomic dislocations on 
the TiO2(H0) surface has a significant effect on its band 
structure. In particular, band edges are sensitive to any 
distortion on the Ti203 construction. For example, sin- 
gle Pt atom adsorption changes the valence and conduc- 
tion band edges and offsets of the reconstructed rutile 
surface when we compare the corresponding DOS struc- 
tures of ar-TiO 2 (110) systems with and without Pt atom 
in Fig. [71 For an isolated Pt, Fermi energy occurs at 0.85 
eV above the VBM due to upper lying gap state. The 
band gap is widened by 0.92 eV in comparison to that of 
the clean reconstructed (1x2) surface. Single Pt atom on 
4x2 cell corresponds to a low concentration coverage. It 
can only partially reoxidize the Ti(iv) defect states, which 
appear in the bottom panel of Fig. [7] for the majority- 
spin component. The states associated to these saturated 
dangling bonds as a result of Pt binding fall into the VB 
at around —2.0 eV. Remaining partially unsaturated de- 
fect states move down to lower energies leading to the gap 
widening. Their DOS exhibits two sharp van Hove sin- 
gularities with dispersions weaker than those of the clean 
reconstructed surface. Non-bonding excess charge leads 
to one and two flat-like filled states just above the VBM 
for majority and minority spin components, respectively. 

A Pt dimer energetically prefers to get adsorbed at the 
trench in interaction with two reactive Ti(iv) cations on 
the added row and with the in-plane Ti5c as shown in 
Fig. [6j). In comparison, the total energy of a 4x2 cell 
with Pt dimer over Ti5c row along [100] is interestingly 
just 0.04 eV higher than the lowest energy geometry. Ev- 
idently, Ti5c is an active site inside the trench. Other 
cases appear to be far less probable. For example, Pt 2 
over and along the Ti 2 03 row is energetically less favor- 
able by 0.92 eV. In Fig.[S]3, the dimer shows an extended 
length of 2.61 A which is considerably larger than its 
isolated value of 2.33 A, essentially due to strong Pt- 
Ti coupling. As a result of the Pt-Ti(iv) interactions 
having 2.44 and 2.56 A bonds, the oxygen over the row 
changes its outward posture and shows an inward align- 
ment. Moreover, Pt at Ti5c site pulls the Ti underneath 
out of its in-plane position by ~0.8 A leading to a bond 
length of 2.28 A. These dislocations distort the Ti 2 03 
group breaking its symmetry along [010]. 

Pt 2 particle brings six majority and five minority flat- 
like impurity states below the Fermi energy as shown in 
Fig. [7] In addition, a slightly more dispersing Pt driven 
state appear just above the VBM for the spin-up com- 
ponent. Fermi energy is determined by the upper lying 
impurity state at 1.17 eV above the VBM. Thus, the band 
gap is larger by 0.75 eV relative to that of ar-TiO 2 (110). 



Evidently, the appearance of partial DOS contribution 
from mainly 5d electrons of the Pt cluster in the vicinity 
of the Fermi energy indicates that Pt 2 has non-bonding 
localized excess charge density which was absent in the 
case of an isolated Pt case. 

We have investigated various probable adsorption 
structures for Pt3 particle. Energetically, Pt trimer fa- 
vors the on-the-row adsorption over an in-trench position 
as shown in Fig. [St- Pt trimer forms in an equilateral tri- 
angle with a 2.78 A and two 2.57 A bonds. Pt3 and Pti 
cases have some common features. Topmost Pt atom 
pulls Ti(iv) atom so that the bond length becomes 2.49 
A similar to that of Pt/ar-TiO 2 (110). The sub-oxygen of 
the row that is denoted by its bond with Ti(iv) as 2.05 
A in Fig. |BJ: is pulled up as in the single Pt case. At 
the other side of the row Pt-Pt bond aligns with [100] 
pushing the two oxygens away so that, that side of the 
row bends outwards. 

The DOS structure of Pt 3 /ar-TiO 2 (110) system dis- 
plays some similarities with that of Pt 2 in that a num- 
ber of impurity states appear below the Fermi energy 
exhibiting dominant Pt 5d character. Evidently, these 
weakly dispersing impurity states are brought by the ex- 
cess charge density localized on the metal cluster. Similar 
to the case Pt 2 , Fermi level is determined by the upper 
lying fiat-like state at 1.27 eV above the VBM that cor- 
responds to a band gap widening of 0.90 eV with respect 
to bare reconstructed surface. On the other hand the 
offsets of the VB edges differ due to adsorbate induced 
different local surface relaxations. Pt3 and a single Pt 
atom stick to the Ti 2 03 group. Therefore the partial 
DOS contribution of bonding 5c? electrons appear in the 
corresponding VBs. Since, the coordination number of 
Pt3 with the added row is relatively larger compared to 
the Pt case, and since Pt3 distorts the row more than 
a single Pt do, bonding Pt3~Ti 2 03 states appear in the 
upper part of the VB. 

The minimum energy structure for the Pt4 case has 
been found to be as shown in Fig|6ji where the trench 
accommodates the cluster. We see that Pt4 adsorption 
builds up on the Pt 2 / ar-Ti0 2 skeleton by additional two 
platinum atoms so that a gas phase Pt4 forms in the 
trench. Inward position of the row oxygen closer to the 
cluster is similar to that of the Pt 2 case. Two adjacent 
vertical Pt-Ti5c bonds become 2.38 and 2.36 A that are 
slightly larger than the value of 2.28 A for the Pt dimer 
on the surface. They reduce the interaction strength by 
relatively better saturating the dangling bonds. Clearly, 
off-planar form of the rhombus gets distorted due to Pt- 
Ti (iv) interactions so that the side bonds are considerably 
extended. As in the case of Pt 2 adsorption the added row 
construction gets pulled towards Pt4 through Pt-Ti (iv) 
interaction. Similar arguments can be made for Pt3 and 
Pt-surface cases. 

In the top panel of Fig. [71 the CB edge is determined by 
empty Pt4 6s states. Eight and seven flat-like occupied 
impurity states appear in the DOS for the spin up and 
down components, respectively, associated mainly with 
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the 5d electrons of Pt4 cluster. The fiat impurity state 
lying 1.59 eV above the VBM sets the Fermi level corre- 
sponding to a band gap larger by 0.28 cV relative to that 
of the bare reconstructed surface. The positioning of the 
cluster-surface bonding states inside the VB is similar to 
Pt2 case which shows common adsorption characteristics. 

Pt adsorbates prefer to bind to Ti atoms on the addcd- 
row model of the reconstructed surface since T12O3 sto- 
ichiomctry leads to oxygen deficiency. Hence, oxygens 
on the added-row become less reactive compared to Ti 
cations. As a trend, odd and even numbered cluster sizes 
show different adsorption characteristics. While clusters 
with odd number of Pt atoms energetically prefer to be on 
the added row terrace, even numbered ones relax into the 
trench at one side of the row. Furthermore, it is interest- 
ing to see that Pt4-surface system has common elements 
with Pt2 case as a building block, although they were 
considered distinctly and their initial structures were not 
relevant in this sense. Pt clusters in the trench breaks 
the symmetrical formation of the added row construc- 
tion along [010] while on-the-row adsorption more or less 
resumes the symmetry. 

IV. CONCLUSIONS 

A systematic analysis of atomic and electronic struc- 
ture of small Pt particles, from monomers to tetramers, 
supported on regular and defect sites of TiO2(110) sur- 
face has been presented based on Hubbard U corrected 
hybrid DFT calculations using Dudarcv's approach. The 
atomic structure of the titania supported Pt clusters de- 
termined by the surface stoichiomctry that constrains the 
charge density. The calculated binding energies of Pt 
particles on the Onishi and Iwasawa model of the (1x2) 
surface happens to be greater than those of the corre- 
sponding adsorption systems on partially reduced (oxy- 
gen vacant) and stoichiometric surfaces due to electron 
derealization from the reduced Ti sites. Interestingly, 
these chemisorbed Pt clusters form geometries similar to 
their gas phase structures except Pt4 at the defect site on 
the partially reduced surface where bridging O-Pt inter- 
action plays a role. In fact, binding energies per Pt atom 



decreases as the Pt cluster size increases due to stronger 
meta-metal coordination. Thus, large Pt particles show 
3D-like nucleation. 

We demonstrated that DFT+U method reproduces 
experimentally observed gap states for the non- 
stoichiometric surfaces where electronic correlation is of 
vital importance. In particular, Ti(iv) 3d states have 
been shown to fall within the band gap of TiC>2(110)- 
(1x2) reconstructed surface. Those were incorrectly pre- 
dicted by pure DFT to be inside the CB giving a metallic 
character, similar to oxygen vacancy states. Interaction 
of small Pt clusters with the TiC>2(110) delocalizes charge 
causing distortion in the geometry around the adsorption 
site. This significantly alters the band edges of the tita- 
nia support and also brings band-gap states depending on 
the cluster size. The position of these defect states show 
strong dependence on the lattice relaxations. Conse- 
quently, adsorbate-substrate coupling leads to significant 
band gap narrowing for bulk terminated and partially re- 
duced surfaces while it gives rise to a gap widening in the 
case of reconstructed surface. No metallization occurs 
for Pt n /Ti02 (n=l-4) systems. These results provide 
good insight into the effect of deposition of small Pt parti- 
cles on the atomic and electronic structure of TiO2(110) 
surfaces with adsorbates. The detailed analysis of the 
differences in DOS upon Pt cluster adsorption puts for- 
ward a sound physical picture for adsorbate-substrate in- 
terface properties. In Pt„-TiO2(110) chemisorption sys- 
tems, the appearance of well localized non-bonding Pt 
impurity states might be useful for catalytic and photo- 
voltaic applications. 
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